The Effect of FOXP3+Regulatory T Cells on Infectious and Inflammatory Diseases

CD4+CD25+FOXP3+regulatory T cells (Tregs) contribute to the maintenance of immune homeostasis and tolerance in the body. The expression levels and functional stability of FOXP3 control the function and plasticity of Tregs. Tregs critically impact infectious diseases, especially by regulating the threshold of immune responses to pathogenic microorganisms. The functional regulatory mechanism and cell-specific surface markers of Tregs in different tissues and inflammatory microenvironments have been investigated in depth, which can provide novel ideas and strategies for immunotherapies targeting infectious diseases.Copyright © 2021. All rights reserved.

Balancing Functional Tradeoffs between Protein Stability and ACE2 Binding in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB Lineages: Dynamics-Based Network Models Reveal Epistatic Effects Modulating Compensatory Dynamic and Energetic Changes.

Evolutionary and functional studies suggested that the emergence of the Omicron variants can be determined by multiple fitness trade-offs including the immune escape, binding affinity for ACE2, conformational plasticity, protein stability and allosteric modulation. In this study, we systematically characterize conformational dynamics, structural stability and binding affinities of the SARS-CoV-2 Spike Omicron complexes with the host receptor ACE2 for BA.2, BA.2.75, XBB.1 and XBB.1.5 variants. We combined multiscale molecular simulations and dynamic analysis of allosteric interactions together with the ensemble-based mutational scanning of the protein residues and network modeling of epistatic interactions. This multifaceted computational study characterized molecular mechanisms and identified energetic hotspots that can mediate the predicted increased stability and the enhanced binding affinity of the BA.2.75 and XBB.1.5 complexes. The results suggested a mechanism driven by the stability hotspots and a spatially localized group of the Omicron binding affinity centers, while allowing for functionally beneficial neutral Omicron mutations in other binding interface positions. A network-based community model for the analysis of epistatic contributions in the Omicron complexes is proposed revealing the key role of the binding hotspots R498 and Y501 in mediating community-based epistatic couplings with other Omicron sites and allowing for compensatory dynamics and binding energetic changes. The results also showed that mutations in the convergent evolutionary hotspot F486 can modulate not only local interactions but also rewire the global network of local communities in this region allowing the F486P mutation to restore both the stability and binding affinity of the XBB.1.5 variant which may explain the growth advantages over the XBB.1 variant. The results of this study are consistent with a broad range of functional studies rationalizing functional roles of the Omicron mutation sites that form a coordinated network of hotspots enabling a balance of multiple fitness tradeoffs and shaping up a complex functional landscape of virus transmissibility.

From Deep Mutational Mapping of Allosteric Protein Landscapes to Deep Learning of Allostery and Hidden Allosteric Sites: Zooming in on "Allosteric Intersection" of Biochemical and Big Data Approaches.

The recent advances in artificial intelligence (AI) and machine learning have driven the design of new expert systems and automated workflows that are able to model complex chemical and biological phenomena. In recent years, machine learning approaches have been developed and actively deployed to facilitate computational and experimental studies of protein dynamics and allosteric mechanisms. In this review, we discuss in detail new developments along two major directions of allosteric research through the lens of data-intensive biochemical approaches and AI-based computational methods. Despite considerable progress in applications of AI methods for protein structure and dynamics studies, the intersection between allosteric regulation, the emerging structural biology technologies and AI approaches remains largely unexplored, calling for the development of AI-augmented integrative structural biology. In this review, we focus on the latest remarkable progress in deep high-throughput mining and comprehensive mapping of allosteric protein landscapes and allosteric regulatory mechanisms as well as on the new developments in AI methods for prediction and characterization of allosteric binding sites on the proteome level. We also discuss new AI-augmented structural biology approaches that expand our knowledge of the universe of protein dynamics and allostery. We conclude with an outlook and highlight the importance of developing an open science infrastructure for machine learning studies of allosteric regulation and validation of computational approaches using integrative studies of allosteric mechanisms. The development of community-accessible tools that uniquely leverage the existing experimental and simulation knowledgebase to enable interrogation of the allosteric functions can provide a much-needed boost to further innovation and integration of experimental and computational technologies empowered by booming AI field.

Stability and expression of SARS-CoV-2 spike-protein mutations.

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants.

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Study of SARS-nCoV2 Indian Isolates Gaining Insights into Mutation Fre-quencies, Protein Stability and Prospective Effect on Pathogenicity

Background: SARS-nCOV-2 is a variant of the known SARS coronavirus family. The mutations in viruses are very rapid and can play a crucial role in the evolution or devolution of the organism. This has a direct impact on "host jumping" and the pathogenicity of the virus. Objective(s): The study aims to understand the frequency of genomic variations that have occurred in the virus affecting the Indian sub-population. The impact of variations translating to proteins and its consequences affecting protein stability and interaction were studied. Method(s): Phylogenetic analysis of the 140 genomes from the India region was performed, followed by SNP and Indel analysis of both CDS and non-CDS regions. This effort was followed by a prediction of mutations occurring in 8 proteins of interest and the impact on protein stability and prospective drug interactions. Result(s): Genomes showed variability in origin, and major branches can be mapped to the 2002 outbreak of SARS. The mutation frequency in CDS regions showed that 241 C >T, 3037 C >T, 2836 C >T, and 6312 C >A occurred in 81.5% of genomes mapping to major genes. Corresponding mutations were mapped to protein sequences. The effect of mutations occurring in spike glycoprotein, RNA dependent RNA polymerase, nsp8, nucleocapsid and 3c protease was also depicted. Conclusion(s): Whilst the mutations in spike glycoprotein showcased an increase in protein stability, the residues undergoing mutations were also a part of drug binding pockets for hydroxychloro-quine. Mutations occurring in other proteins of interest led to a decrease in protein stability. The mutations were also a part of drug binding pockets for Favipiravir, Remdesivir and Dexametha-sone. The work allows analyzing larger datasets to understand mutation patterns globally.Copyright © 2021 Bentham Science Publishers.

SARS-CoV-2 Variants Show a Gradual Declining Pathogenicity and Pro-Inflammatory Cytokine Stimulation, an Increasing Antigenic and Anti-Inflammatory Cytokine Induction, and Rising Structural Protein Instability: A Minimal Number Genome-Based Approach.

Hyper-transmissibility with decreased disease severity is a typical characteristic of the SARS-CoV-2 Omicron variant. To understand this phenomenon, we used various bioinformatics approaches to analyze randomly selected genome sequences (one each) of the Gamma, Delta, and Omicron variants submitted to NCBI from December 15 to 31, 2021. We report that the pathogenicity of SARS-CoV-2 variants decreases in the order of Wuhan > Gamma > Delta > Omicron; however, the antigenic property follows the order of Omicron > Gamma > Wuhan > Delta. The Omicron spike RBD shows lower pathogenicity but higher antigenicity than other variants. The reported decreased disease severity by the Omicron variant may be due to its decreased pro-inflammatory and IL-6 stimulation and increased IFN-γ and IL-4 induction efficacy. The mutations in the N protein are probably associated with this decreased IL-6 induction and human DDX21-mediated increased IL-4 production for Omicron. Due to the mutations, the stability of S, M, N, and E proteins decreases in the order of Omicron > Gamma > Delta > Wuhan. Although a stronger spike RBD-hACE2 binding of Omicron increases its transmissibility, the lowest stability of its spike protein makes spike RBD-hACE2 interaction weak for systemic infection and for causing severe disease. Finally, the highest instability of the Omicron E protein may also be associated with decreased viral maturation and low viral load, leading to less severe disease and faster recovery. Our findings will contribute to the understanding of the dynamics of SARS-CoV-2 variants and the management of emerging variants. This minimal genome-based method may be used for other similar viruses avoiding robust analysis.

In-silico investigation of systematic missense mutations of middle east respiratory coronavirus spike protein.

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe pneumonia-like symptoms and is still pose a significant threat to global public health. A key component in the virulence of MERS-CoV is the Spike (S) protein, which binds with the host membrane receptor dipeptidyl peptidase 4 (DPP4). The goal of the present investigation is to examine the effects of missense mutations in the MERS-CoV S protein on protein stability and binding affinity with DPP4 to provide insight that is useful in developing vaccines to prevent coronavirus infection. We utilized a saturation mutagenesis approach to simulate all possible mutations in the MERS-CoV full-length S, S Receptor Binding Domain (RBD) and DPP4. We found the mutations in MERS-CoV S protein residues, G552, C503, C526, N468, G570, S532, S451, S419, S465, and S435, affect protein stability. We identified key residues, G538, E513, V555, S557, L506, L507, R511, M452, D537, and S454 in the S protein RBD region are important in the binding of MERS-CoV S protein to the DPP4 receptor. We investigated the effects of MERS-CoV S protein viral mutations on protein stability and binding affinity. In addition, we studied all DPP4 mutations and found the functional substitution R336T weakens both DPP4 protein stability and S-DPP4 binding affinity. We compared the S protein structures of MERS-CoV, SARS-CoV, and SARS-CoV-2 viruses and identified the residues like C526, C383, and N468 located in equivalent positions of these viruses have effects on S protein structure. These findings provide further information on how mutations in coronavirus S proteins effect protein function.

Genetic variations affecting ACE2 protein stability in minority populations.

While worldwide efforts for improving COVID-19 vaccines are currently considered a top priority, the role of the genetic variants responsible for virus receptor protein stability is less studied. Angiotensin-converting enzyme-2 is the primary target of the SARS-CoV-1/SARS-CoV-2 spike (S) glycoprotein, enabling entry into the human body. Here, we applied computational saturation mutagenesis approaches to determine the folding energy caused by all possible mutations in ACE2 proteins within ACE2 - SARS-CoV-1-S/ACE2 - SARS-CoV-2-S complexes. We observed ACE2 mutations at residue D350 causing the most stabilizing effects on the protein. In addition, we identified ACE2 genetic variations in African Americans (rs73635825, rs766996587, and rs780574871), Latino Americans (rs924799658), and both groups (rs4646116 and rs138390800) affecting stability in the ACE2 - SARS-CoV-2-S complex. The findings in this study may aid in targeting the design of stable neutralizing peptides for treating minority patients.

Electrostatics in Computational Biophysics and Its Implications for Disease Effects.

This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

The Effects of SARS CoV-2 nsp13 Mutations on the Structure and Stability of Helicase in Chinese Isolates

Objective: Coronavirus Disease 2019 (COVID19) is a viral disease caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS CoV-2). The high mutation propensity of the SARS CoV-2 genome is one of the biggest threats to the long-term validity of treatment options. Helicases are anti-viral targets because of the vital role they play in the viral life cycle. In this study, changes in the protein structure caused by SARS CoV-2 nsp13 mutations were investigated to contribute to the development of effective antiviral drugs. Materials and Methods: Genome data of 298 individuals located in the China location were examined. The mutant model was built using deep learning algorithms. Model quality assessment was done with QMEAN. Protein stability analyses were performed with DynaMut2 and Cutoff Scanning Matrix stability. Changes in substrate affinity were performed with Haddock v2.4. Results: In this study, twenty-eight mutations in nsp13 were identified (23 sense, 5 missense). The changes in protein structure caused by the five missense mutations (Leu14Phe, Arg15Ser, Arg21Ser, Leu235Phe, Ala454Thr) were modeled. The mutations caused a decrease in the stability of SARS CoV-2 helicase (-0.99, -1.66, -1.15, -0.54, and -0.73 for Leu14Phe, Arg15Ser, Arg21Ser, Leu235Phe, Ala454Thr, respectively). The mutations reduced the helicase's affinity to the substrate. The docking scores for wild-type and mutant helicase were -84.4±1.4 kcal.mol-1 and -71.1±6.7 kcal.mol-1, respectively. Conclusion: Helicase mutations caused a decrease in the protein stability and nucleic acid affinity of the SARS CoV-2 helicase. The results provide important data on the development of potential antivirals and the effect of mutation on the functions of viral proteins. © European Journal of Biology.All rights reserved.

IDENTIFICATION OF SURFACE GLYCOPROTEIN MUTATIONS OF SARS-CoV-2 IN ISOLATED STRAINS FROM IRAQ.

Background: The global pandemic of coronavirus disease is a societal, economic, and publichealth crisis that is still underway. The spike glycoprotein of SARS-CoV-2 is one of the primary ingredients for virulence, tissue tropism, and host areas. Aim: This study aimed to determine mutations in the S protein of the Iraqi COVID-19 isolates. Full genome sequences of Iraqi strains were obtained from GISAID. Using statistical saturation mutagenesis and other informatics methods, we investigated 20 sequences of SARS-CoV-2 S protein missense mutation isolates in Iraq selected from NCBI. The following mutations were detected for all the strains under study compared to the wild type: L452R, A522V, E583D and D614G. The number of mutations in the strains was different depending on the location of the state from which the sample was collected The D614G mutation was found in 19 strains. One strain had three mutations, while the other was a wild form strain. The structure of the mutant protein changes dramatically, as does the energy of the atoms concerning the docking position, affecting the protein’s stability. The mutation sites would improve the S protein’s stability. Molecular docking of RBD-ACE2 is affected differently by residues L452R and A522V. (English) [ FROM AUTHOR] Глобальная пандемия коронавирусной инфекции стала длительной кризисной ситуацией для общества, экономики и здравоохранения, которая продолжается и сейчас. Спайк-гликопротеин вируса SARS-CoV-2 является одним из первичных компонентов вирулентности, тканевого тропизма и объектов носительства. Целью работы было определение мутаций S-белка в изолятах от больных COVID-19 в Ираке. Методы: полногеномные последовательности линий вируса в Ираке получали из базы GISAID. Используя статистики сатурационного мутагенеза и другие методики биинформатики, мы изучили 20 последовательностей изолятов SARS-CoV-2 с миссенс-мутацией данного белка, выявленных в Ираке и выбранных из базы данных NCBI. Результаты: во всех линиях вируса, при сравнении с диким типом, были выявлены следующие мутации: L452R, A522V, E583D and D614G. Число мутаций этих линий было различным, в зависимости от места сбора образцов. Мутация D614G была обнаружена в 19 линиях. Одна из линий имела 3 мутации, тогда как другая относились к дикому типу вируса. Структура мутантного белка существенно изменяется из-за энергетических взаимодействий атомов в зоне стыковки, что влияет на стабильность белка. Выводы: стабильность S-белка может изменяться в зависимости от места мутации. Стыковка молекулы RBD-ACE2 нарушается по-разному при заменах аминокислот L452R and A522V (Russian) [ FROM AUTHOR] Copyright of Medical Immunology (1563-0625) is the property of National Electronic-Information Consortium and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full . (Copyright applies to all s.)

The Role of Deubiquitinating Enzymes in Acute Lung Injury and Acute Respiratory Distress Syndrome.

Acute lung injury and acute respiratory distress syndrome (ALI/ARDS) are characterized by an inflammatory response, alveolar edema, and hypoxemia. ARDS occurs most often in the settings of pneumonia, sepsis, aspiration of gastric contents, or severe trauma. The prevalence of ARDS is approximately 10% in patients of intensive care. There is no effective remedy with mortality high at 30-40%. Most functional proteins are dynamic and stringently governed by ubiquitin proteasomal degradation. Protein ubiquitination is reversible, the covalently attached monoubiquitin or polyubiquitin moieties within the targeted protein can be removed by a group of enzymes called deubiquitinating enzymes (DUBs). Deubiquitination plays an important role in the pathobiology of ALI/ARDS as it regulates proteins critical in engagement of the alveolo-capillary barrier and in the inflammatory response. In this review, we provide an overview of how DUBs emerge in pathogen-induced pulmonary inflammation and related aspects in ALI/ARDS. Better understanding of deubiquitination-relatedsignaling may lead to novel therapeutic approaches by targeting specific elements of the deubiquitination pathways.

Biophysical Fitness Landscape of the SARS-CoV-2 Delta Variant Receptor Binding Domain.

Among the five known SARS-CoV-2 variants of concern, Delta is the most virulent leading to severe symptoms and increased mortality among infected people. Our study seeks to examine how the biophysical parameters of the Delta variant correlate to the clinical observations. Receptor binding domain (RBD) is the first point of contact with the human host cells and is the immunodominant form of the spike protein. Delta variant RBD contains two novel mutations L452R and T478K. We examined the effect of single as well as the double mutations on RBD expression in human Expi293 cells, RBD stability using urea and thermal denaturation, and RBD binding to angiotensin converting enzyme 2 (ACE2) receptor and to neutralizing antibodies using isothermal titration calorimetry. Delta variant RBD showed significantly higher expression compared to the wild-type RBD, and the increased expression is due to L452R mutation. Despite their non-conservative nature, none of the mutations significantly affected RBD structure and stability. All mutants showed similar binding affinity to ACE2 and to Class 1 antibodies (CC12.1 and LY-CoV016) as that of the wild-type. Delta double mutant L452R/T478K showed no binding to Class 2 antibodies (P2B-2F6 and LY-CoV555) and a hundred-fold weaker binding to a Class 3 antibody (REGN10987), and the decreased antibody binding is determined by the L452R mutation. These results indicate that the immune escape from neutralizing antibodies, rather than increased receptor binding, is the main biophysical parameter that determined the fitness landscape of the Delta variant RBD.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms.

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Prediction of Residue-specific Contributions to Binding and Thermal Stability Using Yeast Surface Display.

Accurate prediction of residue burial as well as quantitative prediction of residue-specific contributions to protein stability and activity is challenging, especially in the absence of experimental structural information. This is important for prediction and understanding of disease causing mutations, and for protein stabilization and design. Using yeast surface display of a saturation mutagenesis library of the bacterial toxin CcdB, we probe the relationship between ligand binding and expression level of displayed protein, with in vivo solubility in E. coli and in vitro thermal stability. We find that both the stability and solubility correlate well with the total amount of active protein on the yeast cell surface but not with total amount of expressed protein. We coupled FACS and deep sequencing to reconstruct the binding and expression mean fluorescent intensity of each mutant. The reconstructed mean fluorescence intensity (MFIseq) was used to differentiate between buried site, exposed non active-site and exposed active-site positions with high accuracy. The MFIseq was also used as a criterion to identify destabilized as well as stabilized mutants in the library, and to predict the melting temperatures of destabilized mutants. These predictions were experimentally validated and were more accurate than those of various computational predictors. The approach was extended to successfully identify buried and active-site residues in the receptor binding domain of the spike protein of SARS-CoV-2, suggesting it has general applicability.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability.

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Inflammation of sars-cov2-infected human pluripotent stem cell-derived cardiomyocytes and the effects and side effects of remdesivir

Introduction: The global pandemic of the coronavirus 2019 disease (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In addition to respiratory failures, COVID-19 patients exhibited cardiac complications. Studies observed the direct infection and replication of SARS-CoV2 in human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) accompanied by cytopathic effects. However, the underlying mechanisms of SARS-CoV-2-mediated CM death remain poorly understood. In addition, the therapeutic potential of remdesivir (RDV) on CMs has yet to be answered. Methods and Results: We confirmed that SARS-CoV-2 is infectious to and effectively replicates in hPSC-CMs and is cytopathic to hPSC-CMs. We also found that RDV effectively inhibited viral replication at a concentration of 50 nM. RNA-seq analyses demonstrated that expression of immune responsive genes was elevated in SARS-CoV-2 infected hPSC-CMs. Immunostaining and an ELISA assay further revealed formation of inflammasomes and secretion of inflammasome-mediated cytokines, such as IL-1β, IL-18, and IL-6 in SARS-CoV-2 infected hPSC-CMs. RNA-seq analyses showed gene profile changes in SARS-CoV-2 infected hPSC-CMs corroborating with activation of inflammatory signals and cell death pathways. While gene profiles of 0.1 μM RDV-treated SARS-CoV-2-infected hPSC-CMs showed reversal of such changes, a high dose (10 μM) RDV-treated CoV-2-infected hPSC-CMs showed changes in 44% of genes expressed compared to non-RDVtreated CoV2-infected hPSC-CMs. Among those, expression of protein stability related genes, such as genes associated with autophagy and protein ubiquitination increased while expression of antiviral responsive genes decreased. In addition, a high dose of RDV inhibited expression of mitochondrial genes, particularly MitoComplex I and V compositions, which are related to energy production. Conclusions: This study demonstrates that SARS-CoV2 induced inflammasome in hPSC-CMs, which can underlie cardiac damage in addition to direct cytopathic effects. In addition, RDV can reduce inflammasome when introduced early after SARS-CoV2 infection while a high-dose can aggravate cytopathic effects by potential toxicity to mitochondria.

Computational Saturation Mutagenesis of SARS-CoV-1 Spike Glycoprotein: Stability, Binding Affinity, and Comparison With SARS-CoV-2.

Severe Acute respiratory syndrome coronavirus (SARS-CoV-1) attaches to the host cell surface to initiate the interaction between the receptor-binding domain (RBD) of its spike glycoprotein (S) and the human Angiotensin-converting enzyme (hACE2) receptor. SARS-CoV-1 mutates frequently because of its RNA genome, which challenges the antiviral development. Here, we per-formed computational saturation mutagenesis of the S protein of SARS-CoV-1 to identify the residues crucial for its functions. We used the structure-based energy calculations to analyze the effects of the missense mutations on the SARS-CoV-1 S stability and the binding affinity with hACE2. The sequence and structure alignment showed similarities between the S proteins of SARS-CoV-1 and SARS-CoV-2. Interestingly, we found that target mutations of S protein amino acids generate similar effects on their stabilities between SARS-CoV-1 and SARS-CoV-2. For example, G839W of SARS-CoV-1 corresponds to G857W of SARS-CoV-2, which decrease the stability of their S glycoproteins. The viral mutation analysis of the two different SARS-CoV-1 isolates showed that mutations, T487S and L472P, weakened the S-hACE2 binding of the 2003-2004 SARS-CoV-1 isolate. In addition, the mutations of L472P and F360S destabilized the 2003-2004 viral isolate. We further predicted that many mutations on N-linked glycosylation sites would increase the stability of the S glycoprotein. Our results can be of therapeutic importance in the design of antivirals or vaccines against SARS-CoV-1 and SARS-CoV-2.